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Abstract 

In this chapter we treat the quantum damped harmonic oscillator, and study math- 
ematical structure of the model, and construct general solution with any initial condi- 
tion, and give a quantum counterpart in the case of taking coherent state as an initial 
condition. 

This is a simple and good model of Quantum Mechanics with dissipation which 
is important to understand real world, and readers will get a powerful weapon for 
Quantum Physics. 

Keywords : quantum mechanics with dissipation, quantum damped harmonic oscilla- 
tor, general solution, coherent state, quantum counterpart 



1 Introduction 

In this chapter we introduce a toy model of Quantum Mechanics with Dissipation. Quan- 
tum Mechanics with Dissipation plays a crucial role to understand real world. However, 
*E-mail address : fujii@yokohama-cu.ac.jp 



1 



it is not easy to master the theory for undergraduates. The target of this chapter is ea- 
ger undergraduates in the world. Therefore, a good toy model to understand it deeply is 
required. 

The quantum damped harmonic oscillator is just such a one because undergraduates 
must use (master) many fundamental techniques in Quantum Mechanics and Mathematics. 
That is, harmonic oscillator, density operator, Lindblad form, coherent state, squeezed state, 
tensor product, Lie algebra, representation theory, Baker-Campbell-Hausdorff formula, etc. 

They are "jewels" in Quantum Mechanics and Mathematics. If undergraduates master 
this model, they will get a powerful weapon for Quantum Physics. I expect some of them 
will attack many hard problems of Quantum Mechanics with Dissipation. 

The contents of this chapter are based on our two papers [3] and |6J. I will give a clear 
and fruitful explanation to them as much as I can. 



2 Some Preliminaries 

In this section let us make some reviews from Physics and Mathematics within our necessity. 
2.1 From Physics 

First we review the solution of classical damped harmonic oscillator, which is important to 
understand the text. For this topic see any textbook of Mathematical Physics. 
The differential equation is given by 



where x = x(t), x = dx/dt and the mass m is set to 1 for simplicity. In the following we 
treat only the case co > 7/2 (the case u = 7/2 may be interesting). 
Solutions (with complex form) are well-known to be 
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(2.1) 



so the general solution is given by 

x(t) = ax + (t) + ax.(t) = ae-^^^y + ae-^^^y 

= ae -(i^VHF)' + ae'i^V^&y (2.2) 

where a is a complex number. ]f r y/2u is small enough we have an approximate solution 

x(t) « ae _ (a +<w ) t + ae-(?- iw ) t . (2.3) 

Next, we consider the quantum harmonic oscillator. This is well-known in textbooks of 
Quantum Mechanics. As standard textbooks of Quantum Mechanics see [2] and [TT] ([2] is 
particularly interesting). 

For the Hamiltonian 

H = H(q,p) = ^(p 2 + cu 2 q 2 ) (2.4) 

where q = q(t), p = pit), the canonical equation of motion reads 

dH dH 2 

q = — = p, p= -— = -u q. 
op oq 

From these we recover the equation 

q = —uj 2 q <^=^ q + u 2 q = 0. 

See (12.11) with q = x and A = 0. 

Next, we introduce the Poisson bracket. For A = A(q,p), B = B(q,p) it is defined as 

I A B\ - (25) 

c dq dp dp dq 
where {, } c means classical. Then it is easy to see 

{q, q} c = 0, {p,p} c = 0, {q,p} c = 1. (2.6) 



Now, we are in a position to give a quantization condition due to Dirac. Before that we 
prepare some notation from algebra. 
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Square matrices A and B don't commute in general, so we need the commutator 

[A, B]= AB- BA. 

Then Dirac gives an abstract correspondence q — > q, p — > p which satisfies the 
condition 

[<?,<?] = 0, [p,p] = 0, [q,p]=ihl (2.7) 

corresponding to (12. 6p . Here h is the Plank constant, and q and p are both Hermite operators 
on some Fock space (a kind of Hilbert space) given in the latter and 1 is the identity on it. 
Therefore, our quantum Hamiltonian should be 

H = H(q,p) = ^(p 2 + 0J 2 q 2 ) (2.8) 

from (12. 4p . Note that a notation H instead of H is used for simplicity. From now on we 
consider a complex version. From (12. 4p and (12.81) we rewrite like 

H(q,p) = \{p 2 + u 2 q 2 ) = %-{q 2 + \p 2 ) = %-{q - -p){q + -p) 
2 2 ur 2 oj oj 

and 

#(?>P) = ir(9 +-P) = ir W P)(?+-PJ [9,Pj 

= T « P)(? + -P) + - > = < —(q p)(q + -p) + - 

2 I oj oj oj J {2h oj oj 2 

by use of (12.71) . and if we set 

" t = Vi (4 -> a= ^ + i® (2 - 9) 

we have easily 

r + n w rA i _ z' w f 2i.„ .1 oj [ 2% ., 

K a = 75* 9+ -P>9 P = 77T i 9>P ^ = 777- i Xlh 

2n oj oj 2k I oj J 2k I oj 
by use of (12.71) . As a result we obtain a well-known form 

H = u)h(a!a+-), [a,a f ] = l. (2.10) 



Here we used an abbreviation 1/2 in place of (1/2)1 for simplicity. 

If we define an operator N = a) a (which is called the number operator) then it is easy 
to see the relations 

[N,a j ] = a\ [N,a] = -a, [a,a j ] = 1. (2.11) 

For the proof a well-known formula [AB, C] = [A, C]B + A[B, C] is used. Note that aa) = 
a) a + [a, a'] = N + 1. The set {a^, a, iV} is just a generator of Heisenberg algebra and we can 
construct a Fock space based on this. Let us note that a, and N are called the annihilation 
operator, creation one and number one respectively. 

First of all let us define a vacuum |0). This is defined by the equation a|0) = 0. Based 
on this vacuum we construct the n state \n) like 

|n) = ^£|0) (0<n). 



Then we can easily prove 

t|n) = \fn + l\n + 1), a\n) = \/n\n - 1), iV|n) = n|n) (2.12) 
and moreover can prove both the orthogonality condition and the resolution of unity 

oo 

(m\n) = 5 mn , 2j|n)(n| = l. (2.13) 

n=0 

For the proof one can use for example 

a 2 ( a f ) 2 = a(aaV = a(A^ + l)a f = (iV + 2)aa f = (N + 2){N + 1) 

by (12. lip , therefore we have 

(0|aV) 2 |0) = (0|(iV + 2)(N + 1)|0) = 2! (2|2) = 1. 

The proof of the resolution of unity may be not easy for readers (we omit it here). 
As a result we can define a Fock space generated by the generator {a\ a, N} 

{oo oo 
J>n|n> e C°° | ^|c„| 2 <ool. (2.14) 
n=0 n=0 J 



This is just a kind of Hilbert space. On this space the operators (= infinite dimensional 
matrices) a\ a and N are represented as 



/ 



1 

y/2 

\/3 


/ n 



\ 



( 



■ J 





1 

y/2 

\/3 



N = a ] a 



■ J 



by use of fl2T2|) . 

Note We can add a phase to {a, a^} like 



b = e ie a, 6 f 



e a 



t, TV = 6+6 = a f £ 



where 6 1 is constant. Then we have another Heisenberg algebra 



[N,b ] ] = b\ [N,b] = -b, [6,6 f ] = l. 



(2.15) 



Next, we introduce a coherent state which plays a central role in Quantum Optics or 
Quantum Computation. For z G C the coherent state \z) G T is defined by the equation 

a\z) = z\z) and (z\z) = 1. 

The annihilation operator a is not hermitian, so this equation is never trivial. For this state 
the following three equations are equivalent : 

(1) a\z) = z\z) and (z\z) = 1, 

(2) \z) = e zat -~ za \0), (2.16) 
I (3) \z)=e-^EZo^\n). 



The proof is as follows. From (1) to (2) we use a popular formula 

e A Be- A = B + [A,B} + I [A, [A, B}} + • • • 



(A, B : operators) to prove 



e -(«,t-a,) ae *ot-*, = a + z _ 



From (2) to (3) we use the Baker-Campbell-Hausdorff formula (see for example [T7] ) 

e A e B = e A+B+ ^ A ' B]+ * [A ' [A < B]]+ * [B ' [A > B]]+ "\ 
If [A, [A, B}] = = [B, [A, B}] (namely, [A, B] commutes with both A and B) then we have 
e A e B = e A + B +^ A ' B ] = e \[A,B] e A+B ==> . 6 a+b = ^—^[a,b]^a^b (2.17) 

In our case the condition is satisfied because of [a,(v] = 1. Therefore we obtain a (famous) 
decomposition 

The remaining part of the proof is left to readers. 

From the equation (3) in (12.16P we obtain the resolution of unity for coherent states 

r fdxdy^^ = jr\n)(n\ = 1 (z = x + iy). (2.19) 

J J 71 n=0 

The proof is reduced to the following formula 

dxdy 
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^z m z n = n\S mn {z = x + iy). 



The proof is left to readers. See [H] for more general knowledge of coherent states. 

2.2 From Mathematics 

We consider a simple matrix equation 

— X = AXB (2.20) 



X = X(t) 



A 



B 



where 

Xu{t) x 12 {t) 

X2l(t) X 2 2(t) 

A standard form of linear differential equation which we usually treat is 

d 



a ll 0\2 
021 «22 




(It 



x = Cx 



where x = x(t) is a vector and C is a matrix associated to the vector. Therefore, we want 
to rewrite (I2.20p into a standard form. 

For the purpose we introduce the Kronecker product of matrices. For example, it is 
defined as 



A®B 



V 



an Qi2 \ 


"( 




a 21 ^22 J 




021 B 


aii&n on 612 


Onfall 


Ol2&12 


OU&21 011&22 


Ol2&21 


Ol2&22 


a 2 ifeii a 2 ib l2 


022&11 


022&12 


021&21 021&22 


022&21 


022&22 



(2.21) 



/ 



for A and B above. Note that recently we use the tensor product instead of the Kronecker 
product, so we use it in the following. Here, let us list some useful properties of the tensor 
product 

(1) (A x <g> B l ){A 2 <S> B 2 ) = A X A % ® B X B % , 

(2) {A® E){E®B) = A®B = (E®B)(A® E), 



A&E+E&B 



A&EE&B 



(3) e'^" 1 = e"°°~e 

(4) (A ® B) j = A ] ®B ] 



E){E 



e A ® e B , 



(2.22) 



where E is the unit matrix. The proof is left to readers. [9] is recommended as a general 
introduction. 
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Then the equation (I2.20p can be written in terms of components as 



dxn 
dt 

dxi2 
dt 

dx 2 i 
dt 

dx 2 2 
dt 



aii&uxn + aiib 2 \Xi2 + 012^11^21 + «i2&2i^22, 

au6l2Xn + an&22^12 + fll2&12^21 + «12&22^22, 
«21&11^11 + «21&21^12 + «22&11^21 + ^22^21^22, 
a 2 l&12^11 + «21&22^12 + «22&12^21 + Ct22&22^22 



X\\ 



X\2 
X22 



or in a matrix form 



d_ 

dt 



If we set 

X\\ X\2 

X21 X22 

where T is the transpose, then we obtain a standard form 



xn 


\ 


( 


a n b n 


an&2i 


012&H 


G12&21 


X12 






011612 


011&22 


CH2&12 


^12^22 


X21 






a2l6ll 


A21&21 


^22^11 


a 22^21 


X22 


J 


\ 


021^12 


021^22 


&22&12 


022&22 



V 2 / 



X = (x u , X12, X21, ^22)" 



d_ 
~dt 



X = AXB 



d 
dt 



X = (A® B )X 



from 

Similarly we have a standard form 
rf -X = AX + XB 



d 



X = {A®E + E® B T )X 



(2.23) 



(2.24) 



dt dt 

where E T = E for the unit matrix E. 

From these lessons there is no problem to generalize (12.231) and (I2.24p based on 2 x 2 
matrices to ones based on any (square) matrices or operators on J 7 . Namely, we have 



dt 



AXB f t X = (A ® B T )X, 
AX + XB f t X = (A (g) I + I (g) B T )X. 

where / is the identity E (matrices) or 1 (operators). 



dt 



(2.25) 
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3 Quantum Damped Harmonic Oscillator 

In this section we treat the quantum damped harmonic oscillator. As a general introduction 
to this topic see pL] or [T6] . 

3.1 Model 

Before that we introduce the quantum harmonic oscillator. The Schrodinger equation is 
given by 

ih^m)) = HMt)) = C»HN + -)) \*(t)) 

by (12.1 Op (note N = a) a). In the following we use instead of Jr. 

Now we change from a wave-function to a density operator because we want to treat a 
mixed state, which is a well-known technique in Quantum Mechanics or Quantum Optics. 

If we set p(t) = \^f(t))(^f(t)\, then a little algebra gives 

ih—p = [H,p] = [uhN,p] —p = -i[uN,p). (3.1) 

This is called the quantum Liouville equation. With this form we can treat a mixed state 
like 

N 

p=p(*) = X>ii*i(*)><*i(*)i 

3=1 

where Uj > and Ylf=i u j — 1- Note that the general solution of (13.11) is given by 

p(t) = e- iujNt p(0)e iujNt . 

We are in a position to state the equation of quantum damped harmonic oscillator by 
use of ( jSD) . 

Definition The equation is given by 

— p = —i[ua)a, p\ — ^ (a)ap + pa) a — 2apa i ) — — [aa) p + paa) — 2a) pa) (3.2) 
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where /i, v (/i > v > 0) are some real constants depending on the system (for example, a 
damping rate of the cavity mode|l 

Note that the extra term 

— ^ (a)ap + pa) a — 2apa)) — — (aa) p + paa) — 2a) pa) 

is called the Lindblad form (term). Such a term plays an essential role in Decoherence. 



3.2 Method of Solution 



First we solve the Lindblad equation : 







^-p = — — (a)ap + pa) a — 2apa)) [aa) p + paa) — 2a) pa) . 



V 

2 



Interesting enough, we can solve this equation completely. 

Let us rewrite (I3.3P more conveniently using the number operator iV = a)a 

d + + /i + i^,, 7 Ar . ll — v 
—p = papa 1 + ua'pa — (JM p + pM + p) H — p 



(3.3) 



(3.4) 



where we have used aa) — N + 1. 

From here we use the method developed in Section 2.2. For a matrix X = (xij) G M(J-") 
over T 

( 

^00 ^01 ^02 

Xio Xu X\2 

X20 X21 X22 



X 



\ 



we correspond to the vector X E J^dinicJ 7 as 



/ 



X — (Xij) > X — (Xqq, Xqi,Xq2, - - - ! ^lCb ^11) x 12, " ' " ] %20, x 21, x 22i 



(3.5) 



where T means the transpose. The following formulas 



AXB = (A® B T )X, (AX + XB) = (A ® 1 + 1 <g> B T )X 



(3.6) 



lr rhe aim of this chapter is not to drive this equation. In fact, its derivation is not easy for non-experts, 
so see for example the original papers [IS] and [TU], or [T2] as a short review paper 
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hold for A,B,X e M{JF), see flZ25|) . 
Then (13 .4p becomes 



^-p = <j pa g> (a T ) J + z/a T <g> a - (JV ®l + l®iV+l®l) + 1 <g) 1 j p 

^ 2 ^ 1 g> 1 + ua ] g> a ] + pa <g> a - ^-^(iV ®l + l®JV+l<g>l)|p (3.7) 
where we have used a T = from the form (12.151) . so that the solution is formally given by 

p(t) = e ^* e *{^ a+ ® at +Ma®a-^(iVig)l+l(g)Ar+l®l)}^/QN ^3 gN 

In order to use some techniques from Lie algebra we set 

^ 3 = i(jV(g>l + l(g)jV+l® 1), K + = a t ® a f , AT_ = a ® a (iC = ffQ (3.9) 

then we can show the relations 

[K 3 , K+] = K+, [K 3 , KJ\ = -K-, [K + , KJ[ = -2K 3 . 

This is just the su(l, 1) algebra. The proof is very easy and is left to readers. 
The equation (13.81) can be written simply as 

p(t) = e^t^K^K.-^+u)^}^^ (3 1Q) 

so we have only to calculate the term 

e t{vK + +nK-~{n+v)K 3 } ^ (3-11) 

which is of course not simple. Now the disentangling formula in [4] is helpful in calculating 

dm]). 

If we set {k + , k 3 } as 



then it is very easy to check the relations 



[k 3 , k+] = k+, [k 3 , k-] = -k-, [&+, k-] = -2k 3 . 
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That is, {k + , k_, k 3 } are generators of the Lie algebra su(l, 1). Let us show by 577(1, 1) the 
corresponding Lie group, which is a typical noncompact group. 

Since SU(1, 1) is contained in the special linear group SL(2; C), we assume that there 
exists an infinite dimensional unitary representation p : SL(2; C) — U(T ®J-) (group 
homomorphism) satisfying 

dp(k + ) = K+, dp(k.) = K., dp(k 3 ) = K 3 . 

From (13.1 ip some algebra gives 

= p ( e *("fef+A*— by definition) 

= p(e tA ) (3.13) 

and we have 

gtA _ ^t{uk + +pk--(p+u)k3} 



exp < t 



-P e& 



cosh (^t) - sinh (^t) ft sinh (^f) 

sinh (^t) cosh (^t) + ^ sinh (^t) 

The proof is based on the following two facts. 

2 

-p ) v o y v * y \ o i 



and 

oo 



e- = V ^-X" = V A-I 2 " + V — -^-T X 2n+1 (X = tA). 
^ n ^ (2n) ^ 2n + l ! v ; 

n=0 n=0 v y n=0 v y 

Note that 

cosh(x) = = > - — — and sinh(x) = = > — . 

K J 2 ^ 2n ! K J 2 ^ 2n + 1 ! 

n=0 v 1 n=0 v y 
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The remainder is left to readers. 

The Gauss decomposition formula (in SL{2; C)) 



a b 
c d 

gives the decomposition 




(ad — be = 1) 



e tA = 



( ! 7^ sinh (^ f ) 

_^) + i^ sinh (M_^) 



1 cosh(4^ + ^s~inh(^ 

1 



\ 



V 

/ __ 
V 



X 



/ 



i o \ 

cosh(^)+^sinh(^) ) 



x 



and moreover 
/ 



_^_ sin h(ii^t) 



we have 



e tA = exp 



exp 



exp 



\ 



/ 



x 



cosh (4^0 + ^sinh(^ t ) 

y 

' - log (cosh (^) + ^sinh (■ 

V 

/ o o\ 

/ ^i. s inh(^t) 
= eXP ^cosh(^t) + S sinh (^)^ 







log 



X 



exp i 



, log cosh to +! ± 



2\i 



sinh 



eXP l v cosh(^) + S sinh (^) 



(cosh (^t)+^ sinh (^)) 



X 



/ 



\ V 2 I 1 /Li-i/ V 2 / / 

Since p is a group homomorphism (p(XYZ) = p(X)p(Y)p(Z)) and the formula p (e ifc ) = 



14 



X 



^_ sin h(^t) 
G{t) = cosh(^)+^sinh(^) 



x 



e Ldp(k) ^ _ ^ ^ holds we obtain 

/ 4n / -^sinh(^i) 

'<° ) = (cosh (4,7) ^.nh^^V X 

exp ^—2 log ^cosh — ^ + — s i nn ^^~2 — J ) ^P(^3 

As a result we have the disentangling formula 

/ _2£_ sinh (U=Zt) 

eP ^cosh(^t) + ^sinh(^t) " 

exp ^-2 log (^cosh O^f-t) + J^ sinh X 

( ^Mgg „ 

eXP ^cosh(^t) + ^sinh(^t) A - 

by dSH. 

In the following we set for simplicity 

-^-sinh(^t) 
1 ' cosh (tz^t) + ff sinh (^t) ' 



(3.14) 



(3.15) 



Readers should be careful of this "proof, which is a heuristic method. In fact, it is 
incomplete because we have assumed a group homomorphism. In order to complete it we 
want to show a disentangling formula like 



e t{uK++nK—(n+u)K 3 } = e f(t)K +e g(t)Ka e h(t)K- 



with unknowns /(£), g(t), h(t) satisfying /(0) = g(0) = h(0) = 0. For the purpose we set 
A(t) = e *{^++M'--(M+^3} ) B uj _ e f(t)K +e 9(t)K 3e h(t)K-_ 
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For t = we have A(0) = B(0) = identity and 



A(t) = {uK + + fiK. - (/i + u)K 3 }A{t). 



Next, let us calculate B(t). By use of the Leibniz rule 



B(t) 



(fK + )e f ® K +e°® K *e h V K - + e f M K +(gK 3 )e 9 ® K3 e h ® K - + e^ K +e 9 ^ K \hKJ)e h ^ K ~ 
= [fK + + gef K +K 3 e~ fK + + be? K +e° K *K_.e-° K *er* K + } e f® K +e?® Ka e h ® K - 
= [fK + + g(K 3 - fK + ) + he-\K_ - 2fK 3 + f 2 K + )} B(t) 
= {(f-9f + he~ 9 f)K + + (g- 2he' 9 f)K 3 + he~ B K. ) B(t) 

where we have used relations 

e gK 3K e -gK 3 = e -g K _ and e fK+ K_e~ fK+ = K_ - 2fK 3 + f 2 K + . 
The proof is easy. By comparing coefficients of A(t) and B(t) we have 



he 9 = [A, 

g-2he- 9 f = -{n + u) 
f-gf + he~ 9 f 2 = v 

Note that the equation 



he 9 = [A, 
9-2/if = + 

[ f-gf + vf 2 = v 



he 9 = /i, 

9-tyf = + 
[ f + ^ + 1 y)f-^f 2 = u. 



f + ( fI + u)f~^f 2 = u 
is a (famous) Riccati equation. If we can solve the equation then we obtain solutions like 



./ 
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Unfortunately, it is not easy. However there is an ansatz for the solution, G, F and E. That 
is, 

f(t) = G(t), «/(*) = -2 log(F(t)), h(t) = E(t) 
in (I3.15p . To check these equations is left to readers (as a good exercise). From this 

A(0) = 5(0), A(0) = 5(0) ==> A(t) = B(t) for all t 
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and we finally obtain the disentangling formula 

e t{vK + +fiK—bj,+v)K 3 } = e G(t)K +e ~2\og(F(t))K 3e E(t)K. (3.16) 

with (pSJ. 

Therefore ( 13. 8 p becomes 

p(t) = e^r 1 exp (G(t)a f <g> a f ) exp (- log(F(t))(iV ® 1 + 1 ® N + 1 <g> 1)) exp (S(t)a ® a) p(0) 

with ( I3.15p . Some calculation by use of ( I2.22p gives 

pit) = exp (G(t)at <g> a T ) {exp (- \og{F{t))N) ® exp (- log(F(t))iV) T } x 

exp (E(t)a ® (a f ) T ) p(0) (3.17) 

where we have used N T = N and a) = a 1 '. By coming back to matrix form by use of (13.61) 
like 

exp (J57(f)a <8> (atf) p(0) = £ ^j- (a ® (atf) m p(0) 

m=0 

= £ ^r- ( flm « at ) m ) r ) E ^-« m p(o)(« f ) m 

Z — «■ m ! v Z — / m ! 

m=0 m=0 

we finally obtain 

e 2 1 
^ = X 

V^( fl t)» eX p(-log(F(t))iV) J V^ a >(0)(at)-j exp (- \og(F(t))N) 
z — ' n! ' m! 

n=0 L U=0 J 

(3.18) 

This form is very beautiful but complicated ! 



3.3 General Solution 

Last, we treat the full equation (I3.2p 

d is 

—p = —iui(a^ap — pa) a) — — {a)ap + pa) a — 2apa)) — — {aa) p + paa) — 2a) pa) 
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From the lesson in the preceding subsection it is easy to rewrite this as 

^p = ^-iuK + iyK + + f iK_-(fi + iy)K 3 + V^l®l\p (3.19) 
in terms of Kq = iV £g> 1 — 1 £g> iV (note that N T = N) . Then it is easy to see 

[K , K+] = [K , K 3 ] = [K , K-] = (3.20) 
from (I3.9p . which is left to readers. That is, K commutes with all {K + , K 3 , KS\. Therefore 

p(t) = e -i"tKo e t{vK++iiK--(ii+v)K 3 +l^l®l}^Q^ 

= e^r 1 exp (-iutK ) exp (G(t)K+) exp (-2 \og(F(t))K 3 ) exp (E(t)K_) p(0) 
= e^*exp (G(t)K + ) exp ({-iutK - 2 \og(F(t))K 3 }) exp (E(t)KJ) p(0), 

so that the general solution that we are looking for is just given by 



n=0 

oo 



y, E{t^ am f . ^ _ h ,F(t))}N)]a n (3.21) 



mi 

, m=0 



by use of flSTTI) and (gUj) . 
Particularly, if t/ = then 



2 sinh (ft) _ i 



£(t) = ^ — / x = 1 - e 

V ; cosh (ft) + sinh (ft) 

F(t) = cosh (J^tj + sinh (|t) = e*\ 
G{t) = 



from (13.151) . so that we have 



p(t) = e -(f+-M <j 1 , J a m p(0)(a t ) m !> e-(t-- (3.22) 



m! 

, m=0 



from (E2U). 

4 Quantum Counterpart 

In this section we explicitly calculate p(t) for the initial value p(0) given in the following 
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4.1 Case of p(0) = |0)(0| 

Noting a\0) = (<=>■ = (0|a"^), this case is very easy and we have 

n=0 ' \ ' n=0 

because the number operator iV (= a^a) is written as 



P® = f$) J2^r^r\0)(0K = - m J2G(tr\n}(n\ = y^**™" (4-1) 



N = J ^^n\n)(n\ ==>■ N\n) — n\n) 

n=0 

, see for example ( I2.15p . To check the last equality of ( 14. ip is left to readers. Moreover, pit) 
can be written as 

pit) = (1 - Git))e^ G ^ N = e ^-G(t)) e io gG{ t)N^ (4 2) 

see the next subsection. 

4.2 Case of p(0) = \a)(a\ (a G C) 

Remind that \a) is a coherent state given by (I2.16P (a|a) = a\a) (a\a^ = (a\a). First of 
all let us write down the result : 

pit) = e la| 2 e-C^>Ho g G( t ) + io g (i- G(t )) exp |_i ogG ( t ) ( ae -(V-rt")* t + ae'^^a - n) } 

(4.3) 

with Git) in (I3.15p . Here we again meet a term like ( 12. 3 p 



with A = 

Therefore, (14. 3p is just our quantum counterpart of the classical damped harmonic oscil- 
lator. 

The proof is divided into four parts. 
[First Step] From (I3.2ip it is easy to see 

y ^L a m \a){a\^) m = T ^a m \a)(a\a m = V m)l f r \a)(a\ = e E ^ 2 \a)(a\ 
z — ' ml ' ml ^— ' m! 

m=0 m=0 m=0 
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[Second Step] From (I3.2ip we must calculate the term 

e 7Ar |a)(a|e^ = e^e^ 01 " 1 0) <0 1 e~ (afa+ ~" a) 
where 7 = — iujt — log(F(t)) (note 7 7^ —7). It is easy to see 

gCOV — ^ a |Q^ g7-^g aa ^ — l5a g — 7-^g7^|Q^ ^e lN (aa* — aa)e~ lN |g^ ^ae 1 a* — ae ' a |Q^ 

where we have used 

The proof is easy and left to readers. Therefore, by use of Baker-Campbell-Hausdorff 
formula (I2.17P two times 

= e 2 e 2 e (0) = e 2 ^ J |ae 7 ) 

and we obtain 

e^ N \a)(a\e^ N = e^ 2 ^ 7 ^ V 7 )(ae 7 | 
with 7 = — iwt — log(F(t)). 

[Third Step] Under two steps above the equation (13.211) becomes 

e ^ e M 2 (i?( t )-l + e^) - Q(t) n 

n=0 

For simplicity we set z = ae 1 and calculate the term 

G(t) 



Pit) = ^ ^2^-(a^\ae^}(aeV. 



z\a n . 



(« = E^( at )"i*><' 

n=0 

Since \z) = e -l 2 l 2 / 2 e 2at |o) we have 

e- |2|2 ^-=^(a t ) n e zat |0)(0|e 5a a n 

^^2( a t)«| )(o| «l e « 



n=0 



e -|2| 2 e zat 



n! 

. n=0 



e" |2|2 e 2at i ]TG(t) n |n)H J>e 2 

ln=0 

g -|^| 2 e2; at e logG(i)7V e 5a 
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by ( 14. ip . Namely, this form is a kind of disentangling formula, so we want to restore an 
entangling formula. 

For the purpose we use the disentangling formula 

where a, (3, 7 are usual numbers. The proof is given in the fourth step. From this it is easy 
to see 

Therefore (u — > z, v — > logG(t), w — > z) 

- z| 2 (l+log G(t)-G{t)) log G(i) tl lo g G(i)- , r<U\Kr 

(%) = e _|z| e (i-G(t))2 e imH zal +im^ za+losG ^ N 



so by noting 



e —iu)t j 

z = ae 1 = a——- and \z\ 2 = |a| 2 e 7+7 = lal 2 



Fit) F(t) 
we have 

-/ 



a— v. 



is e 2 IrtlSCpm I s ) |a| 2l+ ' og ) G(i) ~ G(t ) ) i°gG(*) Q . e -^t a t I log G(t) ^ mi 1 j C(t)N 
p(f^ — gM l^W - J-Jg 1 1 F(t)^(l-G(t))^ gF(t)(G(i)-l)" e a + F(t)(G(t)-l) Qe a-f-lOg Ur^I^V 



l_L-—h> 



F(t) 

By the way, from (13.151) 

G(t) - 1= -i5y' mm - 1) = ' B(i) - 1= -9r 

and 



1 -G(t)+ log G(t) .-^ J e^T"' 



e-^W — — + log G(t) 



F(t)*(G(f) - l) 2 [ F(t) 

= ^- + e-^logG(t) 

= -(^)-l)+e-^*logG(t) 

we finally obtain 



p(t) = (1 - G(t))e |a| 6 ^ "^^e 



e |a| 2 e-("-" >* log G(t)+log(l-G(*)) g 



logG(i)<ae V at+ae / a-JV 
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[Fourth Step] In last, let us give the proof to the disentangling formula (I4.4p because it is 
not so popular as far as we know. From 



aa) + /3a + 'yN = 70^ a + aa} + (3a 



+ 

V 7/ V 7/ 7 



we have 



6 6 7 6 ^ ^ ' 6 t 



£i£ £ a _s a t 7a t a -at -£ 



Then, careful calculation gives the disentangling formula (I4.4p (N = a) a) 

_«2 £ a _« t 7A r « a t _£ a _M -H| _« t £ a 7 7V s a t _£ a 

e 7eTeT e' e ^ = e t e t e t e -1 e' ei e t 

= e 7 i z g 7 ei e' e ^ 

-(^+^4) - s at 7 7V S-eia 2 t _£ a 

= e 7 7^'eT e'eT e t 

+ e-» -s a t 7A r ^at £ e 7 a -M. a 

— g 7 7 Z 7 Z g 7 g ' g 7 g7 g 7 

_<i a t <* e 7 a t 7A r p£^± a 

— e 7 ^ e t e t e ' e t 
_ e 7 ^ e 7 e ' e t 

by use of some commutation relations 

e sa g tat _ e ^g^ t gsa^ e sa e tN = e tN e seta , e tN e sai = e setat e tN . 
The proof is simple. For example, 

gSa^tat _ ^sagfatg— sa^sa _ ^te aa a) e~ sa ^sa _ ^t(a^+s)^sa _ ^st^ta^sa 

The remainder is left to readers. 

We finished the proof. The formula (14 .3p is both compact and clear-cut and has not been 
known as far as we know. See PQ and [16] for some applications. 
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In last, let us present a challenging problem. A squeezed state \(3) (j3 G C) is defined as 



See for example [I]. For the initial value p(0) = we want to calculate p(t) in ( I3.2ip 

like in the text. However, we cannot sum up it in a compact form like f 14 . 3 [) at the present 
time, so we propose the problem, 
Problem sum up p(t) in a compact form. 

5 Concluding Remarks 

In this chapter we treated the quantum damped harmonic oscillator, and studied mathe- 
matical structure of the model, and constructed general solution with any initial condition, 
and gave a quantum counterpart in the case of taking coherent state as an initial condition. 
It is in my opinion perfect. 

However, readers should pay attention to the fact that this is not a goal but a start- 
ing point. Our real target is to construct general theory of Quantum Mechanics with 
Dissipation. 

In the papers [7] and [S] (see also [13]) we studied a more realistic model entitled "Jaynes- 
Cummings Model with Dissipation" and constructed some approximate solutions under any 
initial condition. In the paper [5] we studied "Superluminal Group Velocity of Neutrinos" 
from the point of view of Quantum Mechanics with Dissipation. 

Unfortunately, there is no space to introduce them. It is a good challenge for readers to 
read them carefully and attack the problems. 
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